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Abstract: A prevailing trend in the stabilization of port-Hamiltonian systems is the assumption 
that the plant and the controller are both passive. In the standard approach of control 
by interconnection based on the generation of Casimir functions, this assumption leads to 
the dissipation obstacle, which essentially means that dissipation is admissible only on the 
coordinates of the closecbloop Hamiltonian that do not require shaping and thus severely 
restricts the scope of applications. In this contribution, we show that we can easily go beyond 
the dissipation obstacle by allowing the controller to have a negative semi-definite resistive 
structure, while guaranteeing stability of both the closed-loop and the controller. 
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1. INTRODUCTION 

In the past decade, passivity-based control (PBC) has 
emerged as a control design method that respects, and 
successfully exploits, the physical structure of a system. 
Using the port-Hamiltonian (pH) formalism, state-space 
control design methods are proposed that lead to con- 
trollers and subsequent closed-loop systems that admit a 
physical interpretation. In this approach the Hamiltonian 
(i.e., the internal stored energy) of the system is the focal 
point of the design method, serving as a Lyapunov function 
> for stability and as a storage function for passivity; see 
(Schaft, 2000) for a basic introduction and (Duindam, 
2009) for a comprehensive summary of the developments 
of the pH framework over the past decade. 

Starting with the so-called energy shaping (ES) routine 
(Ortega, 2001), in which the closed-loop energy is shaped 
using static state feedback, numerous extension have led 
to a variety of control methods that set out to shape the 
system's energy, interconnection and dissipation structure. 
However, only the static full-state feedback methods are 
developed to a form that can be considered generic. Indeed, 
the interconnection and dampingassignment passivity- 
based control (IDA-PBC) method is shown (Ortega, 2008) 
to generate all stabilizing static state feedback controllers 
for pH-systems. 

Existing dynamic output feedback PBC strategies — the 
so-called control-by-interconnection (Cbl) methods — are 
centered around the notion of Casimir functions, which 
statically relate the states of the controller to those of 
the plant (Schaft, 2000; Ortega, 2008). Although these 
Casimir-based control design methods are in some ways 
very attractive, they are not yet developed into a generic 
form. First of all, the input-output structure of both plant 
and controller is assumed to be power-conjugate, i.e., the 
input and output of both plant and controller is assumed 



to constitute a power-port and are thus of equal dimension. 
This can be seen as a drawback, since the input-output 
structure cannot be chosen arbitrarily. The second, and 
more severe drawback is that they are critically hampered 
by the so-called dissipation obstacle (Ortega, 2001). The 
dissipation obstacle dictates that no states can be sta- 
bilized that are subject to pervasive dissipation. Several 
methods have been developed to circumvent the dissipa- 
tion obstacle, see e.g., (Ortega, 2001, 2003; Jeltsema, 2004; 
Ortega, 2008). However, in general, these methods rely on 
changing the actual output of the system and thus require 
a particular input-output structure, which is not always 
possible in practice. 

These considerations call for a more general output feed- 
back method. In this paper, we concentrate on the afore- 
mentioned second drawback and show that by removing 
the passivity constraint on the controller, i.e., by allowing 
for an active controller, the dissipation obstacle can simply 
be resolved while stability of the controller is established 
by using the Casimir relation between the plant and con- 
troller states. The design of output feedback controllers 
for non-collocated input-output systems will be treated 
elsewhere. 

2. CONTROL BY INTERCONNECTION AND 
CASIMIR FUNCTIONS 

In this section, we briefly review the control by intercon- 
nection (Cbl) method applied to port-Hamiltonian (pH) 
plant systems of the form 



y 



[J(0 



R(x)]^-(x)+G(x)u, 



(1) 



where x £ R™, u, y £ M m , J(x) an n x n matrix satisfying 
J{x) = —J T (x), R(x) an n x n matrix satisfying R(x) = 



R T (x), G(x) an n x m matrix, and the Hamiltonian H(x) 
represents the total stored energy. The matrix J(x) is 
usually referred to as the interconnection structure, while 
R(x) captures the resistive structure and is assumed to be 
positive semi-definite. 

If the controller is also a pH system of the form 

d H 

Z=[j c (Z)-R c (Z)]^(0 + G c (Ou c , 



at 



y c = <£(0^(0, 



(2) 



where £ G 



u c ,y c e 



J c (£) an n c x n c matrix 



satisfying J c (£) = — J* (£), R c (£,) an n c x n c matrix satis- 
fying i? c (£) = Rc(Q, G c (£) an n c x m c matrix, and H c (£) 
represents the controller energy, then the interconnection 
of the plant system (1) with (2) via the standard (power- 
preserving) feedback interconnection u = —y c , u c = y 
(assuming m = m c ), yields the closed-loop system 
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(3) 



which is again a pH system. 

The closed-loop system (3) is stabilized if the closed-loop 
Hamiltonian H(x) + H c (£) can be shaped such that it has 
a minimum at the desired equilibrium point x* . The usual 
way to proceed is by restricting the motion of the closed- 
loop system to the subspace 

n = {(x,$)€R n * n °\C(x,Z) = K }, (4) 

with C(x, £) = £ — S(x), where S(x) is assumed to be a 
differentiable function and n G K is some constant, and 
such that the closed- loop Hamiltonian becomes H(x) + 
H c (S(x)+k). This is accomplished if, along the trajectories 
of (3), the functions S(x) are such that 



C(x,0\ 



0. 



The functions C{x, £) are called Casimir functions and are 
independent of the Hamiltonian. We are thus looking for 
solutions S(x) of the partial differential equations (PDE's) 



d T S 
dx 



0, (5) 



~J{x)-R{x) -G(x)Gl{£) 
G c (0G T (x) J c (0 - R c (0\ 

which, under the assumption that R(x) >z and i? c (£) ci 
0, are characterized by the following chain of equalities 
(Schaft, 2000): 

d T S 
dx 



(x)J(x)^(x) 




(6) 




= o, 


(7) 




= 0, 


(8) 


J(x)™(x) 


= -G(x)G%(£). 


(9) 



Unfortunately, the application of the Cbl method is 
severely stymied by the condition (7), which, roughly 
speaking, dictates that the Casimir functions cannot de- 
pend on the coordinates that are subject to dissipation. 
This means that dissipation is admissible only on the 
coordinates of the closed-loop Hamiltonian that do not 
require shaping. For that reason, this condition is referred 
to as the dissipation obstacle (Ortega, 2001). However, 
the dissipation obstacle stems from the assumption that 
both the plant and controller dissipation structures satisfy 
R(x) y and R c {£) h 0. Although these properties 
are necessary to ensure that the plant and the controller 
are both passive systems, they are merely sufficient for 
passivity of the closed-loop system. In fact, the passivity 
assumption of the controller is unduly restrictive, as is 
illustrated in the following example. 

3. A MOTIVATING EXAMPLE 

Consider the RLC circuit shown in Fig. 1. Let <j> denote 
the flux associated to the inductor and q denote the charge 
associated to the capacitor, then the equations of motion 
are given by the pH description 
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where the Hamiltonian H(<fi,q) — ,,, 
the total stored energy. 
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Fig. 1. RLC with pervasive dissipation. 

This circuit was brought forward in (Ortega, 2001) as 
one of the simplest examples of a system that suffers 
from the dissipation obstacle. Indeed, since the equilibrium 
point equals (<p*,q*) = (Lu* /r,Cu*), which is non-zero 
for all u* ^ 0, we need to shape the Hamiltonian in both 
coordinates, but condition (7) dictates that 

= °- (") 

Hence S can only depend on the flux linkage 0, instead of 
both <j) and q as required. 

On the other hand, it is easily shown that the circuit is 
trivially stabilized by a feedforward control of the form 
u = u* . Although this controller is not passive, it docs 



allow for an energy shaping interpretation since by setting 
u = u* the circuit dynamics (10) become 

'4>/L\ , [1 
q/C\ + [o 

<f>/L 
q/C 

<f>/L-u*/r 
q/C-u* 

which can be associated to a shaped Hamiltonian of the 
form 
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such that 
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4. BEYOND THE DISSIPATION OBSTACLE 

In the previous section, it is shown that the dissipation 
obstacle is a direct consequence of the passivity require- 
ment on both the system and the controller. In this section, 
we show that by removing this constraint, the dissipation 
obstacle is easily resolved. The issue that remains then, 
is how the closed loop can be rendered stable. Using the 
Casimir relation between the controller state £ and the 
plant state x, we show that the static state feedback 
interpretation of the closed-loop system solves this issue 
in a straightforward manner. 

Consider again the system of PDE's (5), or equivalently, 
d T S, 



dx 
d T S 
dx 



(x) [J(x) - R(x)} - G c (0G T (x) = 0, (12) 



(x)G(x)G^(0+[M0-Rc(0] =0. 



(13) 



Substituting (12) into the transposed of (13) yields 

— {x) [J(x) - R(x)] -^(x) - [J c (0 + R c (0] = 0, 

which, after separation of the symmetric and skew- 
symmetric part, implies for the relationship between the 
controller and the plant dissipation and interconnection 
structure that 



R c{ = --^ { x)R { x)- { x) 
J c( = ^)J(*)g(*), 



(14) 
(15) 



respectively. It is directly observed that (15) coincides 
with condition (6), whereas (14) coincides with (7) and 
(8) if, and only if, we only allow for dissipation structures 
satisfying R(x) y and i? c (£) h 0. However, if we also 
allow for a non-positive controller dissipation structure, we 
can simply proceed from the closed-loop plant dynamics 
OH . 
dx 

which, by transposing (12), i.e., 

G{x)G T c {H) = -[J{x) + R{x)]—{x), 



x = [J(x) R(x)] ^(x) G(x)G%(0-^(0, 



can be rewritten as 

x=[J(x)-R(x)]—(x) 

no f)J-f 

+ [j( x )+R( x )]—(x)^(0. (16) 
Hence, if J(x) — R(x) is invertible, we can write 
x =[J(x) - R(x)] (j^{x) + [J(x) - R(x)]- 1 

x [J(x) + R(x)] —(x)^(S(x) + «)J , (17) 

which, after suitable choices of the controller Hamiltonian 
H c (S(x) + k), may be interpreted as an energy shaping 
(ES) process such that the closed-loop plant dynamics take 
the form 



x=[J(x)-R(x)]-^(x), 



with 



x) 

x [J(x)+R(x)]—(x)-^(S(x) + k), 



(18) 



^(x) = ^(x)+[J(x)-R(x)}- 1 



satisfying 



(x) = 



d 2 H di 



(x) I , for all x e 



dx 2 \ dx 2 

Furthermore, if R(x) >: 0, and 

dH d d 2 H d 

it follows that x* is a stable equilibrium of (16). 



(19) 



(20) 



In general, stability of (16) does not necessarily imply 
closed-loop stability. However, since S : W l — > W 1 " 
is differcntiable by assumption (as is needed in (5)), it 
is continuous. By virtue of the well-known preservation 
of convergence under continuous mappings (Kolmogorov, 
1999), as x converges to x* , the controller state £ converges 
to £* = S(x*). Hence, stability of x implies stability of £. 
This means that the controller does not need to have a 
positive semi-definite dissipation structure in order to be 
stable in the closed loop. 

Although the method described above extends the tradi- 
tional Cbl method considerably, it is still hampered by the 
assumption that J{x) — R{x) is invertible, and that (19) 
should be satisfied. Both these assumptions are needed for 
a pure ES interpretation of this Cbl method, but are, in 
general, overly restrictive. However, using arguments simi- 
lar to those of the static state feedback IDA-PBC method, 
the Cbl method above can be extended to a dynamic 
output feedback IDA-PBC method. Indeed, starting from 
(16), we then have to look for matrices J d {x) = —J d {x) 
and R d {x) = R d (x), and a function H d : M™ -> M 
satisfying (19), such that 



f) J-f 

R(x)]^(x) 



[J(a 

+ [J(x)+R(x)]^(x)^(S(x) + k) 
[j d (x)-R d (x)]-^(x). 



(21) 



This result effectively generalizes both the standard pas- 
sive ES Casimir-based control method as presented in Sec- 
tion 2, and the extended ES Casimir-based control method 
of (16)— (18), to a dynamic output feedback IDA-PBC 
strategy. Summarizing, we have the following proposition. 

Proposition. Consider the interconnection of the plant (1) 
with the controller (2), assume that (5) holds, and that the 
closed-loop plant dynamics satisfy (21) (resp. (18)). Then, 
if the closed-loop plant Hamiltonian Hd(x) satisfies (20) 
and Rd{x) h (resp. R(x) >; 0), the equilibrium point 
(x*, S(x*)) of the overall system (3) is stable. 



5. A MOTIVATING EXAMPLE (CONT'D) 

Let us return to the RLC circuit of Fig. 1. In Section 3, 
we have seen that the circuit suffers from the dissipation 
obstacle. However, suppose we interconnect, again via 
standard feedback, the circuit (10) with a pH controller 
of the form 

r) M 







'0 -1 " 






1 -1/r 



dH d 



dHo 
8q 



(<M 
(0,9) 



(24) 



Vc 



(22) 



Then, proceeding from (5), we find 
— (4>,q)-G c = 



dS , t . IdS., N „ 

^(M) + -^(M) = o 



dS 



{<f>,q)G c -R c = 



with shaped Hamiltonian 

H d {4, q) = ^ (0 - Lu*/rf + ±(q - Cu*f. 

At this point it is important to emphasize that, although 
the circuit is stabilized by a constant control, the under- 
lying controller dynamics (22) take the form 

i = \{u*?-u*u c , (25) 
y c = -u*. 

Despite the fact that R c < 0, for all u* ^ 0, the controller 
dynamics is asymptotically stable since from (24) we may 
deduce that <j> — > Lu*/r and q — > Cu* , as t — > oo, 
which, since £ = S{(j),q), implies that the controller state 
converges to £* = -u*(q* - <j>* /r). Also note that u c -> 
u* jr implies £ — > 0. 

The Casimir approach thus leads us naturally to the feed- 
forward controller that is known to stabilize the system, 
and in spite of its feedforward character, the controller 
dynamics are instrumental for the construction of the 
controller. 

5.2 Output Feedback Control 



> S((j) q) = G (q — 4>/r) us nex ^ ex Pl°it tne freedom in choosing H c (£) as 



and R c = —-G%, which is non-positive for all r < oo and 
G c ± 0. 

5.1 Recovering the Feedforward Controller 

Now, setting the controller Hamiltonian H c {£) = £, with 
£ = S(cp, q), the dynamics of the controlled circuit becomes 



Hc(Q = r 2 n(* + <ui, 

with £ = S(<f>, q) and dj E R, for i = 1,2, some constants 
to be defined. With this choice, however, it is not possible 
to interpret the control action as an ES process (18) since 



-dq-^ q \ 



1 2/r 
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-1 
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. dq 
dH 
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W>,<z) 



+ 



(-Vc) 



is not a gradient vector field to which we can associate a 
Hamiltonian H d ((f>,q) satisfying (19). On the other hand, 
the closed-loop dynamics of the plant take the form 

'axGcL/r -(1 + ai G c C) 
|_<7J = [ 1 -1/r 
which can be rewritten in pH form as 
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(23) 
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l~dq ^ + ~d~T . 



- ai G c L/r ai G c C/2 
ai G c C/2 1/r 



Rd 



dH n 



dq 



(<M 



(<M 



(26) 



where 



Since ^r{Q = 1, the stabilization problem is trivially 
solved by setting G c = —u* such that (23) is equivalent to 



HM q) = ^ L {4>- La/rf + ^{q~ Caf, 
with a = -r 2 a 2 G c (r 2 + axCC^ 2 - aiG c L)~ x . 



Note that the closed-loop plant dynamics is stable if the 
controller parameters oi and G c are selected such that 
Rd h 0. The desired equilibrium point is determined by 
an appropriate selection of a 2 . 

6. FINAL REMARKS AND OUTLOOK 

In this paper, an extension of the Casimir-based control- 
by-interconnection method is presented that resolves the 
so-called dissipation obstacle. This result is based on 
removing the passivity constraint on the controller and 
using the Casimir relation between the plant and controller 
states to guarantee stability of the closed loop. 

Using a simple RLC circuit example with pervasive dis- 
sipation, two intriguing controller synthesis solutions are 
presented. First, it is shown that this dynamic output- 
feedback based design methodology is able to produce a 
feedforward controller, while the second solution shows 
how the dynamics of the controller are instrumental in 
shaping the energy in the coordinates that are not used 
for feedback. 

Although the results presented in this paper only empha- 
size the possibility of an active controller, future work 
will include conditions under which also active plants are 
allowed, as well as the consideration of pH systems with 
direct feedthrough and output feedback for non-collocated 
input-output channels. 
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